Large-Scale Bulk Motions Complicate the Hubble Diagram 
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We investigate the extent to which correlated distortions of the luminosity distance-redshift re- 
lation due to large-scale bulk flows limit the precision with which cosmological parameters can be 
measured. In particular, peculiar velocities of type la supernovae at low redshifts, z < 0.2, may 
prevent a sufficient calibration of the Hubble diagram necessary to measure the dark energy equation 
of state to better than 10%, and diminish the resolution of the equation of state time-derivative pro- 
jected for planned surveys. We consider similar distortions of the angular-diameter distance, as well 
as the Hubble constant. We show that the measurement of correlations in the large-scale bulk flow 
at low redshifts using these distance indicators may be possible with a cumulative signal-to-noise 
ratio of order 7 in a survey of 300 type la supernovae spread over 20,000 square degrees. 
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Introduction. — The challenge to discover the nature 
of dark energy is pushing all methods and measures 
of cosmology to their limits. The luminosity dis- 
tances to type la supernovae (SNe) which first re- 
vealed the cosmic acceleration 0, 13, 0, are now 
being pursued to obtain tighter constraints on cos- 
mological model parameters 0, H, El Obser- 
vational programs, such as the Supernova Legacy 
Survey (http://www.cfht.hawaii.edu/SNLS/), the Su- 
pernova Factory (http://snfactory.lbl.gov), Essence 
(http://www.ctio.noao.edu~wsne/), the Carnegie Super- 
nova Project (http://cspl.lco.cl/~cspuserl/CSP.html), 
in addition to ongoing efforts by existing groups, are cur- 
rently underway, hoping to achieve ~ 10% constraints on 
the dark energy equation of state parameter. In order to 
decisively advance our understanding, and test for a pos- 
sible time-evolution of the dark energy, a dedicated space- 
based mission is planned as part of the NASA/DOE Joint 
Dark Energy Mission (JDEM). 

The luminosity distance-redshift relation, however, has 
a basic limitation as a tool for cosmology in an inhomo- 
geneous universe. Large scale structures distort the dis- 
tances and redshifts. It is well known that peculiar ve- 
locities of SNe induced by the internal properties of host 
galaxies and clusters contribute a random component to 
distance estimates which can be reduced by averaging 
over many SNe. Furthermore, gravitational lensing of 
SN light reduces the accuracy with which the true lumi- 
nosity distance can be determined for an individual SN 
[cllf^ llollll| . thereby complicating an easy interpretation 
of the Hubble diagram. The effect comes from the slight 
modification of the observed SN flux due to lensing by 
the intervening large-scale structure 0, 0, Q, 0] and 
correlates distance errors of SNe spread over the sky at a 
few degrees or less, due to survey geometries in the form 
of "pencil beams" or long, but narrow strips (lif. Our 
primary concern in this paper is large-scale bulk flows 
|l7| . peculiar motions that are coherent on scales above 
~ 60 Mpc, which correlate individual SN distance esti- 



mates spread over ten or more degrees angular scale. In 
this case the effect comes from the slight Doppler shift- 
ing of both the source and observer, affecting both the 
inferred redshift and the flux, resulting in a non-linear 
correction to the luminosity distance. This correlated 
noise cannot be reduced simply by increasing the sample 
size and is expected to affect the error budget from low 
to intermediate redshifts (z < 0.2). Because the Hubble 
diagram at these low redshifts must be pinned down ac- 
curately in order that we may hope to find a possible time 
variation in the dark energy equation of state E0, it 
follows that accounting for bulk motions is a necessity. 

Fluctuations and anisotropies in luminosity distance 
have been studied previously, with most of the focus on 
formalism |20j, |2l| and the role of gravitational lensing 
(e.g. [22|L The attention has only recently expanded to 
include peculiar motions j^. Our intention is to exam- 
ine the consequences of correlated distortions of luminos- 
ity distances due to bulk motions for the interpretation 
of the Hubble diagram and efforts to extract cosmologi- 
cal information about dark energy. Turning the problem 
around, we will also examine whether low-redshift SNe 
can provide a way to measure large-scale bulk flows. 

Calculational Method. — In order to further quantify 
these statements, we will first summarize the errors in- 
duced by peculiar velocity fluctuations. The effect re- 
sulting from velocities involve two differences: first, the 
inferred redshift is modified by the difference in the veloc- 
ity of the source relative to the observer, projected along 
the line of sight; second, the motion at the observer leads 
to a dipole correction to the distance. In combination, we 
obtain (see, Ref. |20j for details including their equation 
3.15; also [HE!): 
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where n is the unit vector along the line of sight, vsnc 
is the SN velocity, v b s is the velocity of the observer, x 
is the comoving radial distance to the SN, and the prime 
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FIG. 1: Correlations in peculiar velocities, C vv (zi, 22, 0), as a function of two SN redshifts, 21 and 22 and their projected 
angular separation 9. On the left, in panel (a), we consider the correlations as a function of 9 when 2 = z\ — 22 and also with 
Z\ fixed at 0.2 with 22 varied. On the right, in panel (b) we consider the correlations as a function of redshift, with z\ = 22, 
for variety of illustrative 9 values. The horizontal line shows the intrinsic SN measurement error corresponding to 8m = 0.1. 
Peculiar velocities correlate SNe separated at 10 to 100 square degrees on the sky at redshifts around 0.1, but extra covariance 
from SNe at redshifts greater than a few tenths is negligible compared to the intrinsic error. The lensing variance overtakes 
the effect due to correlated motions at redshifts starting at 2 ~ 0.2, depending on the velocity-velocity separation angle. 



denotes the derivative with respect to the conformal time. 
The covariance matrix of errors in luminosity distance is 

Cavy « af nt % + C vv (z h z h %) , (2) 

where af nt is the variance term that affects each distance 
individually {e.g. due to random velocities, or the in- 
trinsic uncertainty in the calibration of SN light curves) . 
C vv (zi, Zj, is defined as the correlation, at redshifts 2, 
and Zj with a projected angular separation of 9^ on the 
sky, due to velocity fluctuations. Eq. J2J) defines the full 
covariance matrix due to peculiar velocities. The covari- 
ance in luminosity distances can be computed, following 
Refs. Hi!, whereby 

27-1-1 2 
C vv (z i ,z j ,9 ij ) = V — — cos Oij—Fi (3) 

£ ' 47T 7T 

even £ 

x J k 2 dkP vv (k, Zi, Zj)ji(k[xi - Xj cos%])jj(/eXj sin%) 

with Fi = (I — l)!!/[2'/ 2 (//2)!]cos7vr/2, and the summa- 
tion is over even values of I. We assume that SNe are 
point sources that trace the linear velocity field, but if 
there is a velocity bias, then the correlations could be 
enhanced. Note that P vv (k 1 Zi, Zj) is the power spectrum 
of velocity fluctuations between redshifts Zi and Zj re- 
spectively, which can be written as 

P vv {k,Zi, Zj ) = D'iz^D'iz^P^i^/k 2 (4) 

where P mm . is the mass fluctuation power spectrum and 
D is the mass growth factor. This form only accounts 



for linear fluctuations at large scales. The variance re- 
lated to velocity fluctuations can be obtained in the limit 
where Zi — Zj and 9ij — > 0. We additionally included 
nonlinear velocities, corresponding to internal motions of 
SNe within halos such as groups and clusters, and found 
that these also do not affect error estimates. This is due 
to the fact that the velocity-induced variance is smaller 
than the intrinsic error, of nt . The effect on the Hubble 
diagram, however, is not negligible since correlations be- 
tween errors are dominated by the large-scale bulk flows 
at low-redshifts. 

In Figure 1(a), we show the luminosity distance covari- 
ance C vv (zi, Zj,9ij) with equal redshifts and also with 
Zi = 0.2 for different values of Zj, as a function of the 
separation angle 0ij. In Figure 1(b) we show the covari- 
ance as a function of z — Zi — Zj for several illustrative 
values of 9ij. For reference, we also plot the variance 
as a function of redshift z and compare it to an intrin- 
sic SN magnitude error of 5m — 0.1, which is the ex- 
pected level to which SN light curves may be calibrated 
in upcoming searches. We note that the recent Supernova 
Legacy Survey (SNLS) has reached an average intrinsic 
error of 0.12 (in magnitudes). Peculiar velocities are a 
concern for SNe separated by angular scales of tens or 
more square degrees as seen in Figure 1(a), and at low 
redshifts, z < 0.2, as seen in Figure 1(b). 

To determine the impact on cosmological parameter 
estimates, we compute the Fisher information matrix 

(5) 

If the errors are uncorrelated in the Hubble dia- 
gram, then the final error on a given cosmological 
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FIG. 2: Expected errors on cosmological parameters due to the peculiar velocity covariances. In Figure 2(a) (left panel) the 
expected errors on a constant dark energy equation of state, w — wo, and the matter density parameter are shown. In Figure 2(b) 
(right panel) the expected errors on the dark energy equation of state parameters, tun and wi, where w = Wo + (1 — a)wi, and 
assuming a prior uncertainty on Qm of 0.01, are shown. We have assumed a survey of 300 SNe out to a redshift of 0.2 over 
10,000 sqr. degrees, such as the Supernova Factory, and another 1700 SNe between redshifts 0.2 and 1.7 in 10 sqr. degrees on 
the sky, such as from SNAP/JDEM. We break the error ellipses to covariances from lensing and peculiar velocities. 



parameter obtained by model fitting is <y- m t/ ' yNjzi)- 
But in the case that there are correlations between 
data points due to bulk flows, the final error is 
close to Omt-\/l + \N(zi) — l]r 2 / \/N(zi) where r = 
Cov(i, j)/ -\/Var(i)Var(j) is the average correlation be- 
tween data points. The limit r — > corresponds to the 
case of uncorrelated errors, but in the limit of perfect 
correlation, r — ► 1, the error remains as dint with no im- 
provement from the number of SNe in the survey. For 
< r 2 < 1, while there is an improvement with increas- 
ing the SNe sample size, in the limit of large numbers, 
the error on an individual parameter will not improve 
beyond rcr int . 

Analysis. — To estimate cosmological parameter mea- 
surement errors, we consider a survey with 2000 SNe, 
similar to the combined Supernova Factory and the 
Supernova Acceleration Probe (SNAP) proposal for a 
JDEM. We distribute 300 SNe uniformly in redshift be- 
tween to 0.2 over an area of 10,000 sqr. degrees, and 
1700 between redshifts 0.2 and 1.7 over 10 sqr. degrees. 
We calculate the covariance matrix of size 2000 by 2000 
obtained by assuming a distribution of separations that 
peaks at roughly one half of the diagonal of the survey 
geometry (assumed to be a square). 

In Figure 2, we summarize our results related to cos- 
mological parameter estimates. Here, we have consid- 
ered the measurement of four parameters, the matter 
density parameter f2 m , the Hubble parameter h which 
can also be considered as an overall normalization to 
the Hubble diagram (and affected by low-redshift bulk 
flows) , and assume a dark energy equation of state given 
by w(a) = wq + (1 — a)w%. Our fiducial test model 
is a cosmological constant plus cold dark matter, with 
Wq = — 1, wi = and matter density il m = 0.3. In 



Figure 2(a), we assume w± = exactly and consider the 
measurement of f2 TO and w — wq. In Figure 2(b), we set 
a prior on Q m with cr(O m ) = 0.01, and consider measure- 
ment of wo and w\ . The error ellipses show the expected 
errors based on which part of the covariance is included. 
The innermost ellipse is the case where covariance is ig- 
nored and only intrinsic noise is included while the out- 
ermost ellipse is the case where both peculiar velocities 
and lensing covariance are taken into account. We can 
see that the velocity correlations dilate the £l m — Wo un- 
certainty by ~ 25% in the case illustrated by Figure 2(a), 
and the wo — w\ uncertainty by ~ 20% in the case illus- 
trated by Figure 2(b). Including the effects of both veloc- 
ities and weak lensing, for which the variance rather than 
covariance between sources is dominant, we see that the 
uncertainties expand by 40%, 150% on wq, f2 m respec- 
tively in case (a), and 50% on wq, w\ in case (b). 

In general, smaller separations at high redshift lead to 
an increase in parameter errors from lensing, while at low 
redshifts correlations at the scales of a few tens sqr. de- 
grees increase the peculiar velocity contribution. A com- 
bination of large area (~ 10,000 sqr. degrees) at z < 0.2 
and a smaller area (~ a few tens sqr. degrees) at higher 
redshifts provides the optimal combination, though co- 
variances are not simply reduced to zero in that case. 

There is some possibility to use the dispersion in the 
Hubble diagram as a measure of peculiar velocity fluc- 
tuations. This is summarized in Figure 3, where we 
plot the angular power spectrum of line-of-sight projected 
velocities as a function of redshift in the form of rms 
fluctuations given by v ims = y/l 2 Ci/2 / Kc. The plotted 
power spectrum is equivalent to the Fourier transform 
of Eq. (3) except that we have not included factors of 
(1 — a/a'x) which relate fluctuations in the velocity field 
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FIG. 3: The rms fluctuations in the line-of-sight projected 
velocities ^l 2 Ci/2nc in km/s as a function of the multipole. 
We show the angular power spectrum at different redshifts. 
The diagonal, dot-dashed line is the expected noise power 
spectrum for velocity measurements for a sample of SNe with 
an intrinsic uncertainty of 0.1 magnitudes and a surface den- 
sity of 50 sr _1 . Along the lines of a survey such as the Su- 
pernova Factory, which expects ~ 300 SNe at low redshifts, 
z — 0.01 — 0.08, and assuming the survey to be spread over 
a fractional sky area / s k y = 0.5, we plot the expected error 
boxes for angular power spectrum measurements binned in 
multipole space, for the mean survey redshift. The velocity 
fluctuations are detected at the cumulative signal-to-noise ra- 
tio ~ 7. 



to that of the luminosity distances estimated with SNe. 
We are assuming that an individual distance estimate, 
combined with redshift and an external estimate of the 
overall normalization of the SN light curve, can be con- 
verted to an estimate of the peculiar velocity [2(|. In 
Figure 3, we also show the uncertainty related to pe- 
culiar velocity fluctuation measurements, cr^/N where 
(T vc i is the intrinsic error in the velocity measurement 
from each SN and N is the surface density of SNe (in 
sr _1 ). Assuming an intrinsic uncertainty of 0.1 mag- 
nitudes, then cr vc i = Gmtcz /2.Y! (in km/s) which at 
z ~ 0.02 is 275 km/s. We also assume no uncertainty in 
the observer's velocity, and that the measurements are 
not limited by uncertainties in cosmological parameters 
such as the Hubble constant or affected by any system- 
atic biases. At low redshifts, surveys such as the Su- 
pernova Factory expect ~ 300 SNe over 2-7T sr so that 
using an estimate of 50 sr -1 for the surface density and 
/sky = 0.5 for the fractional sky coverage, we obtain the 
expected error boxes for binned multipole measurements 
in Fig. 3. The line-of-sight projected velocity anisotropy 
power spectrum is detected with a cumulative signal-to- 
noise ratio of ~ 7 for the noise curve and error bars shown 
in Figure 3. However, this is not a significant detection 
for detailed cosmological parameter estimates. For com- 
parison, unlike the low velocity anisotropy "signal" cap- 
tured by low redshift SNe, it may be possible to study 
clustering statistics of lensing magnification with samples 



of SNe at a redshift beyond unit y w ith signal-to-noise ra- 
tios of order thirty or more [22J, |27J ■ 

While peculiar velocity fluctuations in the Hubble di- 
agram do not provide extra cosmological information, 
there are significant implications for distance estimators 
and cosmological probes. For example, due to the corre- 
lations, a full JVtot x N% t Fisher matrix, as opposed to 
a redshift-binncd smaller version, is required in order to 
obtain cosmological parameter accuracy estimates. Pre- 
sumably this is not a problem since the correct treatment 
of SN calibration uncertainties already requires the full 
iVtot x Mot (or even larger) covariance matrix The 
challenge is significant, however, as one can neither ignore 
small correlations nor assume some arbitrary cosmology 
to estimate covariance among measurements which are 
then used to extract new cosmological parameters. The 
full covariance matrix must be established as a function 
of cosmological parameters to obtain an accurate gauge 
of cosmological parameter uncertainty. 

None of these considerations will deter upcoming 
searches for SNe for cosmological purposes, though a 
careful consideration must be given to account for ve- 
locity fluctuations at low redshifts and lensing effects at 
high redshifts. Since peculiar velocity correlations are 
only significant at z < 0.2, one can potentially ignore low 
redshift SNe when fitting distance data to cosmological 
estimates. In this case, we find that the parameter errors 
are not significantly affected by velocity correlations ex- 
cept that the errors are increased by the fractional factor 
in which the SN sample is reduced. In fact, this increase 
is larger than the case where all SNe are used to estimate 
cosmology, but with a proper accounting of the correla- 
tions. So, instead of simply throwing away data, such as 
low-redshift SNe, it may be best to keep the sample as 
a whole, but develop techniques to account for peculiar 
velocity correlations. 

Since the low, I < 6 multipoles in the velocity 
anisotropy spectrum dominate the SN distance covari- 
ance, if such multipoles can be determined independently 
of the SN measurements then corrections, at least par- 
tially, can be applied to the interpretations of the Hubble 
diagram. If a signal-to-noise ratio of 10 measurement in 
each multipole is adequate for a reasonable correction, 
then independent bulk flow measurements at redshifts 
ranging over 0.01 to 0.1 must involve a source surface 
density of 10 3 sr -1 and an uncertainty in the velocity 
measurement of each object below 100 km s _1 . Such a 
surface density of sources and a velocity error may be 
achievable with cluster studies of the kinetic Sunyaev- 
Zeldovich 29] effect with the upcoming Planck surveyor 
(http://www.rssd.esa.int/index.php?project=PLANCK), 
though foregrounds and internal motions within clusters 
will contaminate bulk flow measurements and reduce 
the overall signal-to- noise ratio levels 30|. Another 
approach will be to consider information from an almost 
all-sky peculiar velocity survey based on low-redshift 
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H' 1 = d/dz[d L /(l + z)} whereby 
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FIG. 4: Correlations in the Hubble constant, C HH (zi, Z2, 9), 
as a function of 6 when Z\ = z%. At redshifts less than 0.05, 
the correlations are generally at the level of few times 10 -4 
and peak at angular scales of 50° , suggesting that the Hubble 
constant should show fluctuations at the level a few percent, 
at most, when divided to patches on the sky at the same angu- 
lar scale. A reliable detection of this few percent fluctuation 
is challenging given the low surface density of SNe expected at 
low redshifts, similar to the detection of velocity anisotropics 
shown in Figure 3. 



galaxy samples. In the past, the IRAS Point Source 
Redshift Catalog (PSCz) has allowed modeling of the 
spherical harmonic moments of the velocity field [3l| 
out to a redshift of 0.02. We encourage the development 
of techniques to use information from such surveys to 
correct the correlations in the low redshift part of the 
Hubble diagram. 

It is interesting to note that the dimensionless fluctu- 
ations of the angular-diameter distance d& due to large- 
scale bulk motions are identical to that for the luminosity 
distance, 



5dA/dA = Sd L /d L 



(6) 



This means our results for velocity covariances apply 
equally to distances based on angular-diameter measure- 
ments. Possible scenarios include the distances to large 
scale structure obtained through baryon acoustic oscil- 
lations [32, distances to galaxies using the Sunyaev- 
Zeldovich effect (for a recent example, see |3j|) or other 
features such as radio lobes [2 UK i an d probes of cos- 
mology using the Alcock-Pazcynski j3(| test which em- 
ploy the angular-diameter distance to a correlation ra- 
dius. The survey details will differ in all cases, so that 
the extent to which velocity correlations of low-redshift, 
wide-separation objects contribute noise will also vary. 

We can also consider dimensionless fluctuations of 
the Hubble constant, inferred from either luminosity or 
angular-diameter distances. Using c?l, we note that 
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The second term on the right includes the correction due 
to the peculiar acceleration. At low redshifts, this is 
equivalent to 5r = —5d L + 5 Z {1 + ^(1 — q)z + 0(z 2 )), 
which we obtain by perturbing the redshift expansion of 
the luminosity distance. Here, 5x is the fractional pertur- 
bation to X = H, dL, zand<5 z = (l+l/z)n-(vsNc~ v bs). 
(See Refs. |2fJ,|23j.) In principle, the deceleration param- 
eter q, also varies on the sky and suffers from correla- 
tions. Similar to fluctuations associated with distance in 
Eq. one can define a covariance for the Hubble con- 
stant anisotropies using the line-of-sight projected cor- 
relation function for the velocity field. This covariance 
is 



C HH (z i ,z j ,0 ij ) = £ 
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where Fi is defined below Eq. (0J. Compared to fluc- 
tuations of the luminosity distance, anisotropies in the 
Hubble constant are larger by a factor of ~ 3 — 5 de- 
pending on the redshift and the deceleration parameter 
(see Figure 4). The increase comes from the correction 
to fluctuations associated with peculiar acceleration in 
Eq. (J2J. As shown in Figure 4, at redshifts between 
0.025 to 0.05, fluctuations in the Hubble constant are 
at most a few percent, given that the correlations are 
C HH {z l ,z J ,9, j ) - O(10~ 4 ) at angular scales of 60 de- 
grees. Detecting such a small fluctuation from a low 
redshift SN survey such as the Supernova Factory, how- 
ever, will be challenging just as velocity fluctuations are 
marginally measurable from SN surveys. This is mostly 
due to the low surface density of SNe expected at low 
redshifts. 

Note that the expressions for the Hubble constant and 
distance fluctuations depend on the line-of-sight source 
and observer velocities separately. Hence, there is the 
possibility of a "Hubble bubble," large fluctuations in 
H, if the motion of the reference frame defined by the 
sources, SNe or large scale structure, does not converge 
to our reference frame, defined relative to the cosmic mi- 
crowave background [23, S EM Ell • For example, a local 
low-density bubble could bias H high by ~ 5% |4lll42l |. 
although observations suggest that the reference frames 
have indeed converged by length scales ~ 50 Mpc/h 
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|43l l44j . Nevertheless, our local motion will induce cor- 
related velocity fluctuations if it is not removed from the 
data properly. 

To conclude, we have investigated the correlated dis- 
tortions of the luminosity distance-redshift relation due 
to large-scale bulk flows and how these correlations limit 
the precision with which cosmological parameters can be 
measured. At low redshifts, peculiar velocities correlate 
errors of type la SNe and prevent a precise calibration 
of the Hubble diagram, relative to the scenario where 
one arbitrarily assumes no correlations so that the er- 
rors decrease by the square-root of the number of SNe. 
The increase in individual error of a distant SN, or the 
correction to variance from the velocity field, however, is 
negligible relative to an expected intrinsic error of 0.1 to 
0.15 magnitudes. These results are consistent with other 
recent calculations on how pecu liar velocities affect cos- 
mological studies with SNe 23]. Turning our arguments 
around, we find that the measurement of large-scale bulk 
flows at low redshifts using SN distance indicators is 
challenging. At high redshifts, weak gravitational lens- 
ing magnification adds an extra dispersion to the light 
curve and increases the individual errors of SN distance 
estimates. For surveys that are concentrated on smaller 
areas on the sky, weak lensing also correlates distance 
estimates, but the increase in individual variances gen- 
erally dominates the error budget. We also considered 
similar distortions of the angular-diameter distance, as 
well as the Hubble constant. 

We thank Caltcch for hospitality, where part of this 
work was completed. R.C. was supported in part by NSF 
AST-0349213 at Dartmouth. 
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